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An approach is proposed enabling to effectively describe, for relativistic heavy-ion collisions, the 
observed deviation from unity of the intercept A (measured value corresponding to zero relative 
momentum p of two registered identical pions or kaons) of the two-particle correlation function 
C(p,K). The approach uses g-deformed oscillators and the related picture of ideal gas of g-bosons. 
In effect, the intercept A is connected with deformation parameter q. For a fixed value of q, the 
model predicts specific dependence of A on pair mean momentum K so that, when |K| J> 500 — 600 
MeV/c for pions or when |K| J> 700 — 800 MeV/c for kaons, the intercept A tends to a constant 
which is less than unity and determined by q. If q is fixed to be the same for pions and kaons, 
the intercepts A,r and Xk essentially differ at small mean momenta K, but tend to be equal at K 
large enough (|K| J> 800 MeV/c), where the effect of resonance decays can be neglected. We argue 
that it is of basic interest to check in the experiments on heavy ion collisions: (i) the exact shape 
of dependence A = A(K), and (ii) whether for |K| > 800 MeV/c the resulting A^ and Xk indeed 
coincide. 



PACS numbers: 25.75.-q, 25.75.Gz, 03.75.-b 



*E-mail: anch@bitp.kiev.ua 
^E-mail: omgavr@bitp.kiev.ua 
"''E-mail: mmtpitp@bitp.kiev.ua 



1 



Hadron matter under intense conditions of high tem- 
peratures and densities is extensively studied in relativis- 
tic heavy-ion collisions (RHIC). Models and approaches 
used to describe the processes in the reaction region 
are examined by comparing predictions of these mod- 
els with experimental data on single-, two- and many- 
particle momentum spectra. Two-particle correlations 
are known to carry information about the space-time 
structure and dynamics of the emitting source Usu- 
ally, the study of correlations occurred in RHIC assumes 
that (i) particles are emitted independently (completely 
chaotic source), and (ii) finite multiplicity corrections can 
be neglected. Then, correlations reflect the effects from 
(anti)symmetrization of the amplitude to detect identical 
particles with certain momenta, and the effects due to fi- 
nal state interactions (FSI) of detected particles between 
themselves and with the source. The FSI, known to de- 
pend on the structure of emitting source, as well jgrovide 
certain information about source dynamics, see 

The correlation function is defined as || 



C(k a ,k b ) = 



P 2 (k a ,k b ) 
P (k Q ) Pi (k b ) 



(1) 



Pi (k) and P2(k a ,k b ) being single- and two-particle 
probabilities to registrate particles with definite mo- 
menta. 

For identical particles, the two-particle wave function 
appears to be a symmetrized (antisymmetrized) sum of 
single-particle states (chaoticity assumption), namely 

i>7 a 7b(Xa,Xb,t) = -j=[lp la (x a ,t) ^ (x&, t) + 

+ e M ^ 7a (x b ,i)V 7b (x a ,i)]- (2) 

Here a = (a — n) for identical bosons (fermions). The 
indices 7 a ,7 b of the 1 -particle wave functions label com- 
plete sets of 1 -particle quantum numbers. Below, we refer 
our consideration of two-particle correlations to identical 
bosons (pions, kaons, etc.). 

In the absence of FSI, for chaotic source, the correla- 
tion function can be expressed as follows 



C(p,K) = 



(3) 
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with the 4-momenta K and p defined as K — h(k a + 
kb) , p = k a — kb . The source function S(X, K) (single- 
particle Wigner density) is defined by the emitted single- 
particle states ip 7 (x) at freeze-out times || (for details of 
the derivation see Appendix in j^j), namely 



S(X,K) = 
= [ d 4 xe lK - x 



(4) 



Using hermiticity of the density matrix p 77 ' one easily 
shows that S(X,K) is real. Due to this, at zero relative 
momentum one has C(0, K) = 1 + cosa. From the latter 
relation it is seen that the boson correlation function at 
a = should approach the exact value 2 as the relative 
momentum — ► 0. But, as observed from the very first 
experiments and up to most recent data, the measured 
correlation function never attains this value at p = 0. 
To remove this discrepancy, the correlation function of 
identical bosons is usually represented in the form 



C(p,K) = l + Xf( P ,K) 



(5) 



with A drawn from a fit to the data, A = 0.4 - 0.9, 
f(p, K) commonly taken as Gaussian, so that /(p = 
0, K) = 1. The deviation of A from unity in RHIC can be 
interpreted by production of secondary pions from reso- 
nance decays going outside the fireball. Presence of long- 
lived resonances results in an increase of measured source 
size and life-times HQ . 

Trying to explain experimental data, with Eq. (^|) in 
mind (admitting that a, besides or tt, can as well 
take other values) it is natural to relate the parameter A 
with the effective angle a and to get the reduction factor 
by means of cos a. The correlation function (^), which 
is measurable quantity, possesses the obvious symmetry 
a — > —a, hence there is no contradiction in taking the 
wave function in the form (^|). One can argue that the 
two-particle wave function of a boson pair released from 
a dense and hot environment effectively acquires an ad- 
ditional phase. The drawn phenomenon can be ascribed 
to properties of the medium formed in RHIC, which thus 
exhibits a non-standard QFT behavior through the con- 
sidered correlation functions. Adopting that correlation 
function approaches 1 + A when the relative momentum 
p — > 0, we construct an effective model capable to mimic 
real physical picture. For this, we use ]Bj q-deformed com- 
mutation relations and the techniques of q-boson statis- 
tics which result in a partial supression of the quantum 
statistical effects in many-particle systems. 

As expected, usage of appropriate g-algebra allows to 
reduce the treatment of complex system of interacting 
particles to that of a system of non- interacting ones, at 
the price of more complicated (deformed) commutation 
relations. The deformation parameter q is viewed as 
an effective (not universal) parameter which efficiently 
encapsulates most essential features of complicated dy- 
namics of the system under study. For example, in the 
application of g-deformed algebras to description || of 
rotational spectra of superdeformed nuclei, q is also non- 
universal and assumes a particular value for each nu- 
cleus. In the context of hadron theory the quantum 
(q-deformed) algebras also proved to be useful. Such a 
usage significantly improves description of hadron char- 
acteristics both in hadron scattering |9]~|Tl|] (nonlinearity 
of Regge trajectories) and in the sector of hadron masses 
and mass sum rules Q. 

In this letter we propose to employ, for the system of 
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pions or kaons produced in RHIC, the ideal g-Bose gas 
picture based on so-called g-bosons. Physical meaning or 
explanation of the origin of ^-deformation in the consid- 
ered phenomenon depends on whether the deformation 
parameter q is real or takes complex value, say, a pure 
phase factor. Here we do not deal with the diversity of 
all possible reasons (not completely unrelated) for the 
appearance of q-deformed statistics. Let us mention only 
that the compositeness of particles (pions, kaons) may 
also lead |13[ to a (/-deformed structure with real q. 



We use the set of g-oscillators defined as 14 1 



-Sijbj 



[Nub}] 



[Ni,Nj]=0 



[bi,b j } = [bl,b\} = , 



b i b]-q-^b)b i = S ij q N ^ 



(6) 



g-distribution function (10) turns out to be real, owing 
to dependence on q through the sum q + q . 

Note that the g-distribution (9) already appeared in 
refs. [fl6| (moreover, a two-parameter generalization of 
distribution (10) based on appropriate two-parameter de- 
formed version of the relations (6) is also possible [p"8|). 



f (k) 



I - Bose-Einstein 

II - q-Bose-Einstein 

III - Maxwell-Bolt zmann 
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In the Fock type representation, with the "g-bracket'' 
defined as [r] q = (q r — q~ r )/(q — <7 _1 ), we have 



b\h = [Ni_ 



(7) 



The equality b\bi — Ni holds only if q = 1 (no-deforma- 
tion limit). It is required that either q is real or 



q = exp(i0) 



< 9 < 7T 



(8) 



Below, just this exponential form will be adopted for q 
(compare with the phase a in Eqs. (0),(3)). 

For a multi-pion or multi-kaon system, we consider the 
model of ideal gas of g-bosons (IQBG) taking the Hamil- 
tonian in the form fllq,[l7| 



(9) 



m" + kf, and 



with i labeling energy eigenvalues, u>i = \J 
Ni defined as above. This is the unique truly noninter- 
acting Hamiltonian with additive spectrum. Also, we as- 
sume discrete 3-momenta of particles, i.e. the considered 
system is contained in a large box of volume ~ L 3 . 

Basic statistical properties, as usual, are obtained by 
evaluating thermal averages (A) = Sp(Ap)/Sp(p) , p = 
e~@ H , now with the Hamiltonian (9). Here /3 = 1/T and 
the Boltzmann constant is set equal to 1. 

With b\bi = [N t ] q and q + q~ l = [2} q = 2cos6», the 
q-deformed distribution function is obtained as 



(bth) = 



z 2 ^ - 2 cos(0)e0"" + 1 



(10) 



If q — > 1 it yields the Planck-Bose-Einstein distribution, 
as should, since at q = 1 we return to the standard system 
of bosonic commutation relations. Although the defor- 
mation parameter q is chosen in the complex form (8), the 



FIG. 1. The (/-distribution function versus momentum 
(curve II), in comparison with the quantum Bose-Einstein 
(curve I) and Maxwell-Boltzmann (curve III) distributions. 
The inputs are: T = 120 MeV, m — m^. Curve II corresponds 
to the deformation angle 8 = 24° . 

The shape of the function f(k) = (b^b)(k) from (10) 
corresponding to the gas of pions modelled by IQBG is 
given in Fig. 1. As seen, the q-deformed distribution 
function lies completely in between the well-known two 
curves (standard Bose-Einstein distribution function and 
classical Maxwell-Boltzmann one), demonstrating that 
the deviation of q-distribution (10) from the quantum 
Bose-Einstein distribution goes in the "right direction" 
towards the classical Maxwell-Boltzmann one, which re- 
flects a decreasing of the manifestation of quantum sta- 
tistical effects. 

Analogous curves for q-distribution functions can be 
given with other fixed data. For kaons, because of their 
larger mass and higher empirical value A ~ 0.88 of the 
intercept (which corresponds to a smaller deformation in 
our model), such a curve should lie closer (than pion's 
one) to that of the Bose-Einstein distribution. 

Now let us go over to our main subject of two-particle 
correlations. Calculation of the two-particle distribution 
corresponding to the g-oscillators @ yields 



{b\b\hh) = 



2cos6» 



2cos(26>)e 



1 



(11) 



Then, the desired formula for the intercept A.; = + 1 = 
{b\b\bibi) I '{(b\bi)) 2 of two-particle correlations follows, 

K + l= 2 ™ 9 t +1 -™ 9 ) 2 (12) 



£f + 2(l-cos 2 9)U 

with ti = cosh(Pui) - 1. Note that (11), (12) are 
real functions since, like (10), these depend on the q- 
parameter (8) through the combination \(q + q^ 1 ) = 
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cos 9. Below, we drop the subscript i and ignore the 
above assumed discreteness of momenta. 

The quantity A can be directly confronted with em- 
pirical data. In the limit q — > 1, from Eq. (12) the in- 
tercept Abe = 1 proper for Bose-Einstein statistics is 
reproduced. This corresponds to Bose-Einstein distribu- 
tion, see (10) at 9 = 0. On the other hand, at 9 = n/2, 
Eq. (12) coincides with the value Afd = — 1 proper for 
Fermi-Dirac statistics. The two cases, Abe and Afd, are 
seen in Fig. 2 as the only two points where all the dif- 
ferent curves (the continuum parametrized by w — (3uj) 
merge and the dependence on momentum and temper- 
ature disappears. From the continuum of curves, there 
exists a unique limiting asymptotic one A = 2 cos 9 (or 
A = — 1 + 2 cos 6*), which corresponds to w — > oo, i.e. to 
low temperature or to large momentum. 

Let us discuss implications of Eq. (12). Solving it for 
cos# at a fixed A = Ai, we obtain the deformation angle 
as the function 9 = 0(Xi, K, T, m). Fig. 2 illustrates the 
properties of intercept (correlation strength) A treated 
from the standpoint of g-deformation, i.e., on the base of 
Eq. (12). Note that the continuum of curves A = A(cos#) 
parametrized by w = 0u> divides into three classes ( "sub- 
continua") given by the intervals: (i) < w < Wo, (ii) 
Wq < w < w' 0l and (hi) w' < w < oo. The two "critical" 
values wo = wb — 0.481 and w' = wd — 0.696 fix 
the curves B and D respectively. The curves A, C, E are 
typical representatives of the classes (i), (ii), (iii). All 
the curves from classes (i), (ii) possess two extrema; the 
curve D is unique since its extrema degenerate, coincid- 
ing with the point of inflection. This enables us to define 
the range of "small deformations" 7 sma u for the variable 
9: from 9 = (no deformation) to the value yielding 
minimal A, A m j n ~ 0.33, implied by the "critical" value 
w' Q = wd- On the interval I sma ii the intercept A mono- 
tonically decreases with increase of 9 (or 1 — cos 9, the 
strength of deformation) . 
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FIG. 2. The (1 — cos #)-dependence of the intercept A 
of two-particle correlation, given by Eq. (12). The curves 
A,B,C,D,E correspond to the values wa = 0.3, wb = 0.481, 
wc = 0.58, Wd = 0.696, we = 2.0 of the dimensionless vari- 
able w = f3ui. 

All the type (ii)-(iii) curves at 9 ^ lie below the line 
A = 2 — largest possible correlation attainable in Bose- 



Einstein case (the curve B contains, besides 9 = 0, single 
point at certain value of 1 — cos 6*, where A = 2 is also 
attained). For the type (ii) curves we restrict ourselves 
with Ismail, ignoring moderately (and very) large defor- 
mations. The class (i) contains "irregular" curves for 
each of which there exist g-deformations yielding corre- 
lation strengths exceeding A = 2. Consider special values 
of physical variables T, |K|, which provide the peculiar 
values Wo = wb — 0.481 and w' Q = wd — 0.696 (recall 
that w = v / m 2 + K 2 /T). With m^' ) = 139.57 MeV 
and lowest mean momentum of the pion pair |K| = 0, we 
get the respective two lower bounds for the temperature: 
T = 290.0 MeV and Tq = 200.5 MeV. Compare these 
values with that for a typical curve from class (iii): at 
w = we = 2.0 (the curve E) for pions with |K| = wc 
get the lowest temperature Te = 69.8 MeV. 

The results similar to presented above (starting with 
Eq. (6)) were also given |5j for g-oscillators with defining 
relation aa^ — qa^a = 1 introduced by Arik and Coon ||. 

Theoretical approaches to RHIC are aimed to find 
an adequate description for the non-equilibrium state 
formed during the collision, and the g-boson technique 
enables us to treat the non-stationary hot and dense mat- 
ter effectively as a "noninteracting ideal gas" . To deter- 
mine the g-parameter corresponding to actual state of 
the hot medium, we propose a way of extracting a useful 
information from the two-particle correlations. In a more 
general context, by this we propose and develop an effec- 
tive picture for the two-particle correlations in RHIC. 

Let us specially emphasize a remarkable feature exhib- 
ited within our model: 

At w — > oo (i.e. for very low temperature at fixed mo- 
menta or very large momenta at fixed temperature) , see 
Fig. 2 and relevant comments, we come to the equality 

A = 2cos6»-l (T->0 or |K| -> oo) (13) 

for the employed q-oscillators (6). This means a direct 
link A <-> q whose explicit form is A = [2] q = 2cos#. 

Finite temperature and momenta, on the other hand, 
become non-trivially involved in the relation between A 
and the parameter q as seen from Eq. (12) and Fig. 2. 

In Fig. 3, we present the dependence, implied by 
Eq. (12), of intercept A on the momentum |K| exem- 
plifying this for pions at T = 120 MeV (and at T = 180 
MeV) with four curves corresponding to the fixed val- 
ues 4°, 9°, 15°, 24° of the deformation angle 9. Anal- 
ogous dependence for the case of kaons is presented in 
Fig. 4. Each curve has its own asymptote given by 
Eq. (13) and lying beneath the line A = 1. Recall that, if 
the unique cause for reducing of the intercept would be 
decays of resonances (commonly accepted idea), all the 
curves would tend to the value A = 1 in their asymp- 
totics. In this respect, on the contrary, our model clearly 
suggests possible existence of some other reasons, than 
just the decays of resonances, for the deviation of A from 
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unity. This property probably reflects complicated dy- 
namics of strongly interacting system prior to freeze-out. 
The particular kind of the behaviour of A with respect to 
|K| (implied for a fixed deformation), is a direct conse- 
quence of our model. 
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FIG. 3. The intercept A versus pion momentum |K| for 
fixed 9: A) 9 = 4°, B) 9 = 9°, C) 9 = 15°, D) 9 = 24°. The 
inputs are T = 120 MeV - solid curve, T = 180 MeV - dashed 
curve, and m = m n . 
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FIG. 4. The intercept A versus kaon momentum |K| for 
fixed 9, with the same as in Fig. 3 values of the parameters 
9, T, and now m = tuk- 

Let us summarize main points of the present letter. 
In contrast to the commonly accepted opinion that the 
intercept A of the two-pion correlations should asymp- 
totically approach the value one, we predict attaining by 
A the value less than one at large pion (kaon) pair mean 
momentum |K|. This effect can be verified in the current 
experiments which take place on the SPS at CERN and 
on the Relativistic Heavy Ion Collider at BNL. There- 
fore we propose to pay special attention, in these ex- 
periments, to the behaviour of the intercept A at large 
enough pion (kaon) pair mean momentum, that is, |K| 
in the range up to 500 — 600 MeV/c for pions (up to 
700-800 MeV/c for kaons) in the fireball frame. If our 
prediction is confirmed, this means that A encapsulates 
certain memory effects: it carries important information 



about the fireball dynamics, especially concerning proba- 
ble quark-gluon plasma phase, hadronization, other non- 
stationary processes which took place in the fireball be- 
fore freeze-out. 

As follows from our model, at large enough pair mean 
momenta of the two pions or two kaons, the dependence 
of intercept A on temperature as well as on the particle 
mass (and energy) is washed out. Due to this, A T and Xk 
should practically coincide already at |K| ~ 800 MeV/c, 
although at much smaller mean momenta the intercepts 
A„- and Xk significantly differ from one another (compare 
the corresponding curves in Fig. 3 and Fig. 4, labelled by 
the same 9). Let us recall that this remarkable feature 
is the result of asymptotical reducing of the full T-, K- 
, and m-dependence contained in Eq. (12) to the direct 
link between the intercept A and the deformation param- 
eter q (with q = expi0), given by Eq. (13). Therefore, 
we put forward two principal points for experimental ver- 
ification: 

1) The shape of dependence of intercept A on the pair 
mean momentum as well as its asymptotics (constant 
value of At,-, Xk less than unity at large |K|); 

2) Coincidence of the intercepts for different pairs of par- 
ticles at large enough mean momentum, for instance, the 
regime with A^ ~ Xk should manifest itself starting al- 
ready from |K| ~ 800 MeV/c. 

Verification of these issues is of fundamental interest. 
Indeed, the first point will give answer to the question: is 
the q-dcformcd description (q-boson field theory) really 
capable to mimic the unusual properties of "hot" and 
"dense" hadron matter? In particular, can we adopt this 
description as an efficient tool to deal with noncquilib- 
rium, short-lived system of pions or kaons at high tem- 
perature and density? The second point of verification 
will clarify whether one can regard the deformation pa- 
rameter q as a kind of quantity, which may be referred 
to as one of basic characteristics of the hot, short-lived 
system (besides "temperature" and chemical potential). 

We hope that the answers to these intriguing questions 
will be deduced in the nearest future in the experiments 
which take place in CERN and BNL. 
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